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Abstract 



The effects of a possible rotation of the galactic dark halo on the calculation 
of the direct detection rates for particle dark matter are analyzed, with special 
attention to the extraction of the upper limits on the WIMP-nucleon scalar 
cross section from the experimental data. We employ a model of dark halo 
rotation which describes the maximal possible effects. For WIMP masses 
above 50 GeV, the upper limit exclusion plot is modified by less than a factor 
of two when rotation is included. For lighter masses the effect can be stronger, 
suggesting the necessity to develop specific models of halo rotation in order 
to provide more accurate conclusions. 
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I. INTRODUCTION 



The possibility to detect Weakly Interacting Massive Particles (WIMPs) distributed in 
the halo of our Galaxy has been a major issue in the last years, since these particles could 
provide the amount of dark matter necessary to explain many observed dynamical properties 
of galaxies, clusters and of the Universe itself. Different kinds of possible signals have been 
identified and looked for, in order to outline the presence of WIMPs in our Galaxy. These 
signals are usually referred to as "direct" and "indirect" detection rates. Direct detection 
refers to the possibility to measure a WIMP-nucleus interaction in a low-background de- 
tector, while indirect detection relies on the measurement of WIMP annihilation products: 
photons, antiprotons and positrons produced by the annihilation in the galactic halo, or 
neutrinos coming out of the Earth or the Sun where WIMPs may have been accumulated 
as a consequence of gravitational capture. It is remarkable that the present sensitivity of 
the different experiments is already at the level of the predicted rates for specific WIMP 
candidates, like the neutralino, which represents one of the most interesting and studied 
cold relic particles |J. 

The calculation of the different detection rates depends not only on the particle physics 
properties of the WIMPs interactions, but also on the characteristics of the galactic halo 
where the WIMPs are distributed. Direct detection rates and upgoing-muon fluxes at neu- 
trino telescopes, which both rely on the WIMP elastic scattering off nuclei, depend on the 
WIMP matter density p and velocity distribution /©(f) at the Earth position r in the 
Galaxy. In particular, the dependence of the signals on p is a linear one. The other indirect 
signals (photon, antiproton and positron fluxes) have a stronger dependence on the matter 
distribution, since they are proportional to the square of the matter distribution function 
(DF) p(r) integrated over the effective region of production and propagation of the anni- 
hilation products. On the contrary, this kind of signals are essentially independent on the 
details of the velocity DF, since the annihilating WIMPs are almost at rest and corrections 
due to their velocity dispersion are negligible. 

Detailed estimates of the detection rates would require specific and accurate models of 
the galactic halo able to provide a reliable WIMPs DF g(x, v) (not necessarily separable in 
phase space, i.e. g(r,v) = p{f)f{v)). Unfortunately, detailed halo models are not available 
at present, mainly because the constraints obtained from astrophysical observations are 
not stringent enough to restrict different possibilities. The most important observational 
constraint is provided by the flatness of the rotation curves at large radii. Although the 
available data on our Galaxy do not provide a compelling evidence of a flat rotation curve, 
this feature is observed in a large number of spiral galaxies and therefore it looks reasonable 
to assume its validity also for our Galaxy. 

The standard and simplest model of the dark galactic halo, which is compatible with a 
flat rotation curve, is the so-called isothermal sphere. This model relies on the two basic 
assumptions of spherical symmetry and thermal equilibrium, which find a strong support 
in the argument of "violent relaxation" introduced by Lynden-Bell 30 years ago 0. In this 
model, the DF is separable into a matter density distribution p(r), which has a r~ 2 behaviour 
at large radii, and into a Maxwell-Boltzmann (MB) velocity DF f(v) ||. Although such a 
model gives a divergent total mass and therefore an appropriate cut-off has to be introduced 
at large radii, its range of validity has been tested at least in the inner parts of many 
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galactic systems. Moreover, since it represents a simple and reasonable approximation, in 
the absence of a more detailed model it is widely adopted to describe the dark halo of our 
Galaxy. However, many different models are known to be consistent with flat rotational 
curves. For instance, models which describe non-spherically symmetric or flattened halo 
distributions have been discussed ||. In these models, the specific form of p(f) differs 
from the standard isothermal sphere matter DF, especially at small radii, entailing quite 
large uncertainties on the local value p & . A comprehensive numerical study which takes 
into account a large number of models indicates that the local value of the non-baryonic 
dark matter density falls in the (rather conservative) range 0.1 < p < 0.7 GeV cm -3 
0]. Contrary to the matter DF, the specific form of the velocity DF f(v) has been much 
less investigated. Modifications to the standard MB velocity DF are known HH, but the 
problem of determining the correct form of the distribution of the WIMP velocities in the 
halo has no clear and simple solution at present, both theoretically and observationally. 
The velocity DF is required to be consistent with a given p(r) but this, in general, does not 
determine f(v ) in a unique way. 

The calculation of the WIMP detection rates is usually performed by using the stan- 
dard isothermal sphere model. However, modifications in the isothermal model can affect 
the detection rates, introducing uncertainties in the theoretical predictions and in the ex- 
traction of the experimental limits on the WIMPs parameters. The effects induced on the 
detection rates by a modification in the matter DF are simple to take into account, since 
the dependence of the detection rates on p(r) can be factorized. Specifically, the physical 
range of p & quoted above implies an uncertainty of about a factor of 7 in the evaluation of 
the direct detection rates and in the neutrino fluxes [|IJ (it has to be remarked that this large 
factor reflects a rather conservative attitude). Even larger uncertainties affect the indirect 
rates from WIMP annihilation in the halo, since in this case a modification in the matter 
density profile can strongly affect the integral of p 2 (r) over the effective production region 
of the signal Contrary to the case of the matter DF, a modification of the standard 

MB velocity DF would affect the direct detection rates and the indirect rates at neutrino 
telescopes in a much more involved way. This is because the dependence of these rates 
on f(v) is through a convolution of f(v) with the differential WIMP-nucleus cross section. 
Since the WIMP-nucleus scattering depends on the relative velocity of the WIMPs with 
respect to the detector nuclei, a potentially significant effect could be due to a bulk rotation 
of the halo. This would necessarily modify the WIMP phase-space DF with respect to the 
standard MB form. 

In this paper we wish to discuss the possible effects induced by a halo rotation on the 
direct detection rates, with special attention to the ensuing consequences on the determi- 
nation of the upper limits on the WIMP-nucleus cross section from the experimental data. 
A calculation of the direct detection rates in the case of a rotating halo has been addressed 
in Ref. ||, where it has been concluded that the maximal effect of rotation leads to a 30% 
effect on the total detection rates for a Ge nucleus, in the case of an ideal detector with no 
threshold. However, when considering a real detector the behaviour of the differential rates 
at threshold and the detector characteristics are crucial in determining the experimental 
limits on the WIMP-nucleus cross section . Therefore, we explicitly take into account the 
features of running detectors, such as thresholds, quenching factors and energy resolution, 
in order to estimate the largest uncertainties induced by a possible halo rotation in a con- 
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fident way. To this aim, following Ref. [[| we model the galactic rotation as described by 
Lynden-Bell in Ref. [0, where the maximally rotating velocity DF compatible with a given 
mass distribution has been derived, on the ground of purely kinematical arguments. Even 
if Lynden-Bell's model of halo rotation may not represent a situation which is realized in a 
physical halo, we consider it useful to bracket the size of the effect of halo rotation on the 
direct detection rates. 

The plan of our paper is the following. In Sect. II we briefly describe the calculation of 
the direct detection rates in the presence of halo rotation. In Sect. Ill we discuss our results 
for Ge, Nal and Xe detectors, taking into account the most recent experimental data of 
the different Collaborations. Finally, in Sect. IV we draw our conclusions. An Appendix is 
added, where we report the analytical expressions of the relevant part of the direct detection 
rates which contain the details of the velocity DF in the case of the standard non-rotating, 
maximally co-rotating and maximally counter-rotating haloes. 



II. DIRECT DETECTION RATES 

The interaction of a WIMP of mass m x with a detector produces the recoil of a nucleus 
with energy E R of the order of few to tens keV. The recoil energy can be measured by means 
of various experimental techniques with different nuclear species. At present, experiments are 
running with Ge, Nal, Xe, CaF 2 , Te02 detectors [pTH 20| and other nuclei are currently under 



investigation. The relevant quantity to be calculated and compared with the experimental 
measurements is the differential detection rate 



dE 



lW N T P± [ dvf(v)v^(v,E R ) (1) 
m v J dE R 



where Nt is the number of the target nuclei per unit of mass, p x is the local WIMP matter 
density, v and f(v) denote the WIMP velocity and velocity DF in the Earth frame (v — \v\) 
and da/dE R is the WIMP-nucleus differential cross section. The nuclear recoil energy is 
given by E R = m 2 ed v 2 (1 — cos 8*)/rriN, where 9* is the scattering angle in the WIMP-nucleus 
center-of-mass frame, m N is the nuclear mass and m red is the WIMP-nucleus reduced mass. 
Eq.(|l]) refers to the situation of a monoatomic detector, like the Ge detectors. For more 
general situations, like for instance the case of Nal, the generalization is straightforward. 
The differential WIMP-nucleus cross section can be expressed as 

^ F\q) (2) 



dE R ££ ax 

where do is the point-like total WIMP-nucleus cross section, E R ax is the maximum value 
of E R and E(q) denotes the nuclear form factor, expressed as a function of the momentum 
transfer q 2 =| q \ 2 = 2rriNE R . 

The nuclear form factor depends sensitively on the nature of the effective interaction 
involved in the WIMP-nucleus scattering. To be definite, in the following we will consider 
the case of a WIMP-nucleus scalar interaction, since this is the one which is currently 
accessible to the present sensitivity of running detectors flll,|21"| . In this case we use the Helm 
parameterization of the scalar form factor 
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F(g) = 3 exp --s g 3 

V 2 / 

where s ~ 1 fm is the thickness parameter for the nucleus surface, r = (r 2 — 5s 2 ) 1 / 2 , 
r = 1.2 A 1 / 3 fm, A is the nuclear mass number and ji(qro) is the spherical Bessel function 
of index 1. 

In the isothermal halo model the velocity DF is a Maxwell-Boltzmann distribution in 
the galactic rest frame. Taking into account a finite escape velocity, its expression is the 
following: 



where the normalization factor is 

N = 



2 

erf(z) =zexp(— z 

7T 



2\ 



(5) 



In the previous Eqs., z 2 = 3t> 2 sc /(2w 2 ) and w denotes the root mean square velocity. In the 
isothermal sphere model, w is related to the asymptotic value v ^ of the rotational velocities 



by the simple relation w = J f Uqo- The measured rotational velocity of the Local System at 



the Earth's position is v(r@) = 220 ± 20 km sec 1 [23| and remains almost flat (to roughly 



15%) between 4 Kpc and 18 Kpc. Identifying this value with one gets the estimate 
w = 270 ± 25 km sec" 1 . 

In order to evaluate the WIMP-nucleus interaction-rate, Eq. @ has to be transformed to 
the rest frame of the Earth, which moves through the Galaxy with a velocity v Q = 232 ± 20 
km sec -1 in the azimuthal direction (this value for v & takes into account the motion of the 
solar system with respect to the Local System). Therefore, the velocity v of the WIMPs, 
as seen in the Earth's frame, is related to their velocity in the Galactic frame i^ 3,1 by the 
following set of transformation equations: v$ = v^ 1 — v & in the azimuthal direction, and 
v± = v 5_ al and vr = fff 1 in the vertical (_L) and in the radial (R) direction along the galactic 
plane. 

By means of the previous definitions, the differential rate can be written in the form 



dR _ AT _ p x m N a p2<>J2> 

2 

^red 



-N T ^^F^)l(v min} v e} v esc ) (6) 



where the function I(v min , %, f esc ) contains all the details of the integration of the velocity 
DF f(v) in the Earth's frame: 



l(v min ,v Q ,v csc ) = f dv^-= f vdvf(v). (7) 

J V -'Wmin(B R ) 

In Eq.(0) we have defined: 

f(v) = 2tt / dcos9 f(v,cos9) (8) 

J (cos 9) m \„ 
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where (cos#) m i n and (cos#) max depend on v, v & and v esc . In Eq.(^) v m i n (Eji) = 
[m N E R / '(2m^ ed )) 1//2 . Moreover, we have explicitly considered that particles which possess 
velocities greater than the escape velocity f esc are not bounded to the halo. The value of 



v esc is somewhat uncertain: w esc = 600 ± 200 Km sec" W% - A low value of v esc can sizably 
affect the detection rates, especially at low WIMP masses. In the next Section, we will use, 
as a reference value, v e sc = 650 Km sec -1 . 

Eq.(^) describes a galactic halo which does not possess a bulk rotation. In order to 
analyze the effect of a possible rotation of the isothermal sphere, we consider a class of 
models discussed by Lynden-Bell [Hi] which, for any given mass distribution, describe the 
fastest rotating steady state by means of the following recipe: 

U(v^) = I /( " SaI) < > ° (9) 

o vT l < o 1 ' 



'-^ = l<o (10) 

The LB model corresponds to an un-relaxed system with the maximal rotation compat- 
ible with a given mass distribution. There is no indication, neither theoretical nor obser- 
vational, that this model can be realized in physical galactic haloes. Nevertheless, since it 
provides the largest rotation effect, we choose to use it in order to estimate the maximal 
modification induced by galactic rotation to the direct detection rates and to the extraction 
of the upper limits on the WIMP-nucleus cross section. 

The analytic expressions for X(t> min , t> , t> esc ) are given in the Appendix for the non- 
rotating model of Eq.(|j) and for the models of Eq.([|) (maximal co-rotation) and Eq.flTUD 
(maximal counter-rotation). The three DFs f{v), for the non-rotating, co-rotating and 
counter-rotating cases, are plotted in Fig.l. In the co-rotating situation the WIMPs have a 
bulk rotation in the same azimuthal direction as the Earth. Therefore, the relative velocity 
between the WIMPs and the detector is, on average, reduced. On the contrary, for a counter- 
rotating halo the relative average velocity increases and, moreover, there is a lower velocity 
cut-off corresponding to v Q . The reason why all the three curves in Fig.l cross at the same 
point is a feature of the particular choice of the distribution functions in Eqs. (|9|)-(|I0D and 
does not reflect any general property of rotating models. 

Eq.(|(J) represents the differential rate for an ideal detector. In order to compare the 
calculated rates with the measured ones, we have to express Eq.(|6|) as a function of the 
electron-equivalent energy E ee (which is actually measured, instead of E R ) |E5J. The quan- 



tity E ee is simply proportional to the nuclear recoil energy through the quenching factor 
Q, i.e. E ee = QEr. Moreover, the energy resolution of the detector has to be taken into 
account. This is obtained by means of the convolution 

dR K f°° ( (E ee -E) 2 \ dR 



dE ee V2^r(E ee ) Jo * V 2r2 ( E oc) J dE 

where K is the normalization factor and is given by K = 2/[l + erf (E cc / \^2r(E cc ))]. The 
resolution function r(E) can be expressed as a function of the energy as 
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(12) 



where the energy E is expressed in keV. The constant parameters a and b depend on the 
detector and are determined experimentally. The values of the quenching factors Q and of 
the resolution parameters a and b are reported in Table 1 for the detectors considered in 
Sect. III. The threshold energies E^ of the same detectors are also given in Table 1. 



III. RESULTS 

We start our analysis by discussing the Ge detectors. Fig. 2 shows the differential rate 
on a Ge detector as a function of the electron-equivalent energy E ee . The solid lines refer 
to WIMPs of masses m x = 60 GeV (upper solid line) and m x = 20 GeV (lower solid line). 
For the upper line we used, as a reference value, <jq = 10~ 9 nbarn and for the lower ones 
o"o = 10~ 10 nbarn. The other parameters used in the calculation of the rates are: v Q = 232 
km sec -1 , t> esc = 650 Km sec" 1 and p x = 0.5 GeV cm -3 . The differential rates for a 
maximally co-rotating halo are plotted as dashed lines (the upper one refers to m x = 60 
GeV, the lower one to m x = 20 GeV). The results for a counter-rotating halo are reported 
as dotted lines. 

For low energies, the co-rotating model gives a larger rate as compared to the non- 
rotating model, since low values of E ee mainly correspond to low WIMP velocities in the 
local frame, where the co-rotating DF is enhanced with respect to the non-rotating model, 
as is shown in Fig.l. On the contrary, for higher values of E ee , the co-rotating rate becomes 
smaller than the non-rotating one, and the counter-rotation situation gives the highest rates. 
The difference among the three curves increases with the energy, the one corresponding 
to a co-rotating halo rapidly diverging from the other two. We also notice that, for a 
fixed value of E ee , the difference between the rotating and non-rotating situations is more 
pronounced for lighter WIMPs. This property has a direct influence on the determination of 
the exclusion plot in the WIMP-nucleus cross section <jq vs. m x plane, since the exclusion 
plot is obtained by comparing the calculated rate with the 90% C.L. upper limit on the 
counts of the detector. Therefore, because of a possible halo rotation, we expect a larger 
uncertainty for lighter WIMPs in the determination of the exclusion plots. 

In Fig. 2 the 90% C.L. upper limits of two representative Ge experiments are also plotted: 
the solid histogram refers to the Neuchatel experiment ||12||, the dot-dashed histogram is 



obtained from the Twin experiment [|Tj]. Because of the relatively fast decrease of the 
differential rates as a function of E ee , the most stringent limits on the WIMP parameters 
are usually obtained from the energy bins closest to threshold energy. For the Neuchatel 
experiment, the energy threshold of the detector is 1.5 keV. In the case of Twin, the detector 
energy threshold is 4 keV, but the most stringent limits are provided by the counting rates 
in the energy bins above the Gallium peak (which is clearly visible as an increase in the 
counting rate around E ee ~ 10 keV). This fact gives an effective threshold for WIMP searches 
of about 12 keV, denoted by a vertical dashed line, similar to the Ge/Heidelberg experiment 
T3| . Fig. 2 shows that for Ge experiments with a low energy threshold the effects of halo 



rotation are less important than for the case where the threshold energy is high. This is 
especially true for a counter-rotating halo and high WIMP masses. 
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We are now in the position of determining the upper limit on the WIMP-nucleus cross 
section as a function of m x by comparing the calculated and the experimental rates. Actually, 
the direct detection technique measures the product (p x x cr ) between the WIMP local 
density p x and the WIMP-nucleus cross section cr . The single parameters p x and a cannot 
be disentangled in a direct detection measurement. Therefore, we report our upper bounds 
in terms of the product £er, where the WIMP local density is parametrized as a fraction £ of 
the local total dark matter density p & (£ < 1) and we choose as a reference value p & = 0.5 
GeV cm -3 . As it was discussed in Ref. |], in the case of scalar coupling it is possible to report 
the results as upper limits on the WIMP-nucleon cross section (instead of a WIMP-nucleus 
one), which is more suitable for the comparison among different experiments, especially 
when they make use of different target nuclei. The WIMP-nucleon scalar cross-section is 
defined as |L] 

(nucioon) _ / 1 + m x /mN \ CT 
ascalar " {l + m x /mj A* 

where m p is the proton mass. 

Fig. 3 shows the 90% C.L. upper limit on the quantity £cr s ^ 1 ^ 1 r eon ' ) as a function of the 
WIMP mass, for the case of a non-rotating halo (solid line), of maximal co-rotation (dashed 
line) and maximal counter-rotation (dotted line). The plot is the convolution of the most 
stringent limits of all the presently available results from Ge detectors |T2|^T5f . In the WIMP 
mass range reported in Fig.3, the upper limit for masses below 25 GeV (for non-rotating 
and counter-rotating halo models, 40 GeV for co-rotation) is provided by the Neuchatel 
experiment and for higher masses the limit comes from Twin. Fig.3 shows that the effect 
of rotation of the halo can strongly affect the low mass region in the co-rotating situation. 
The reason for this behaviour has been previously discussed in relation with Fig. 2, where 
it was shown that the detection rates for lighter WIMPs are more affected by a possible 
halo rotation. On the contrary, for masses above 50 GeV the effect is contained below a 
factor of two, when the co-rotation case is compared to the non-rotating halo result. The 
difference between the two exclusion plots is smaller for heavier WIMPs. As reference values, 
the difference is of the order of 50% for masses around 100 GeV and is reduced below 20% 
for m x > 300 GeV. For counter-rotating models, the modification of the exclusion plot 
with respect to the non-rotating case is always relatively small, not exceeding a factor of 
two in the whole mass range. For WIMP masses above 100 GeV, the magnitude of the 
effect is very similar to the case of co-rotation. We have to remind at this point that the 
effect considered here refers to situations of maximal rotation of the galactic halo. Plausible 
models of rotating haloes lie somewhere in between the two maximal cases discussed here, 
probably much closer to the non-rotating case than to these extreme situations. Therefore, 
our results have to be considered as maximal possible effects of galactic rotation on the 
determination of the exclusion plots. 

Let us now discuss the case of Nal detectors. The large mass low-background Nal 
detector of the DAMA/Nal Collaboration currently provides the most stringent upper limit 
on the WIMP-nucleon scalar cross section ]16| (except for a narrow window around m x ~ 15 
GeV, where the Ge detectors are more sensitive Q.) Fig. 4 shows the differential rate for 
Nal detector as a function of the electron equivalent energy E ee . The solid line is the rate 



8 



calculated for m x = 60 GeV and cfq = 10~ 9 nbarn in the case of a non rotating halo. 
The dashed line refers to the co-rotating case and the dotted line to the counter-rotating 
situation. The histogram is the 90% CL. upper limit from the DAMA/Nal Collaboration 
16fl . The ensuing exclusion plot is shown in Fig. 5. Also for the Nal detector, the co-rotating 



case deviates significantly from the non-rotating situation for relatively low masses. On the 
contrary, for counter-rotating haloes the exclusion plot remains close to the result for a 
static halo. The relative deviation R = [^(T^^ ^/[^^^"^jjuja ro t (i stands for co- or 
counter- rotation) between the extracted upper limits in the case of rotation with respect to 
the non-rotating case, are shown in Fig.6. The dashed (dotted) line refers to a co-rotating 
(counter-rotating) halo. We observe that also in this case, for WIMP masses larger that 50 
GeV the effect of rotation affects the exclusion plot by less than a factor of 2. This effect is 
actually smaller for heavier WIMPs: for m x > 80 GeV the exclusion plot uncertainty is of 
the order of 20-30%. 

Finally, we show in Fig. 7 the effects of halo rotation on the exclusion plot for a 129 Xe 
detector. Due to the high nuclear mass, Xe detectors are in principle more sensitive to 
higher WIMP masses than the Nal and Ge ones. Moreover, the quenching factor of Xe 
detectors is larger than the one of Ge and Nal detectors, and this also shifts the sensitivity 
of Xe detectors to higher WIMP masses. The DAMA/Xe Collaboration recently reported 
the results of the analysis on an improved statistics of 1763.2 Kg x day obtained with an 
enriched liquid Xe scintillator flSfl . The exclusion plot obtained from the data of Ref. [18[ 



is plotted in Fig. 7, for the conservative value of Q=0.44. Similar to the cases previously 
discussed, for a counter-rotating halo the deviation is always smaller than a factor of two. 
The same situation happens for co-rotating models and WIMP masses larger than about 50 
GeV. Regardless of the model of rotation, for m x > 200 GeV the deviation is always smaller 
than 25%. 



IV. CONCLUSIONS 

In this paper we have investigated the effect induced by a possible rotation of the galactic 
halo on the rates of WIMP direct detection. In particular, we have discussed the implication 
of halo rotation on the determination of the exclusion plots on the WIMP-nucleon cross 
section for different detectors, namely Ge, Nal and Xe ones. The rotation of the halo 



has been described by using a model |10[ which corresponds to a situation where the halo 
possesses the maximal rotation compatible with a given mass DF, which for simplicity we 
have chosen to be that of the isothermal sphere. 

We found that the exclusion plots obtained from the data are affected by less than a 
factor of 2 in the case of counter-rotating models. The same size of uncertainty occurs 
also for the co-rotating models, when the WIMP mass is larger than about 50 GeV. For 
lighter WIMPs and co-rotation, the exclusion plots are modified by a larger amount. We 
have to remind that, due to the particular model of halo rotation which we have employed 
here, these are expected to be maximal effects. For specific physical rotation models, the 
effect of halo rotation will be plausibly smaller. We can therefore conclude that, at least for 
WIMP masses greater than about 50 GeV, the determination of the exclusion plots from 
the experimental data are affected by an uncertainty smaller than a factor of two due to 
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the possibility that the galactic halo rotates, independently on the specific model of halo 
rotation. We notice that recent preliminary data from accelerators indicate that the lower 
limit on the mass of the most plausible WIMP candidate, the neutralino, is m x ~ 30 GeV for 
low value of the susy parameter tan/3, and m x ~ 45 GeV for tan/? > 3 |27|. Therefore, for 
this dark matter candidate, the uncertainty on the exclusion plot due to a possible rotation 
of the halo is expected to be relatively small. The situation is different for lighter WIMPs. 
In this case, it would be required to develop specific models of halo rotation in order to 
obtain more accurate conclusions. 
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APPENDIX A: 

In this Appendix we report the analytical expression of the function X(v min , v®, v esc ) 
which enters in the calculation of the differential rate for direct detection Eq.(|6]), for a 
velocity DF which is Maxwellian in the galactic frame (see Eq.@). The cases of a non- 
rotating, maximally co-rotating and maximally counter-rotating halo models are given. In 
all the following expressions we make use of the dimensionless variables: x 2 = 3v 2 /(2w 2 ), 
T] 2 = 3vq/(2w 2 ), z 2 = 3v 2 sc /(2w 2 ), x min = (3mE R /Am 2 eA w 2 ) 1/2 , where v Q denotes the Earth 
velocity in the galactic frame, t> esc is the escape velocity, w is the root mean square velocity 
of the Maxwellian distribution in the galactic rest frame. In the definition of the three 
distributions, the same value of w has to be used, since they correspond to the same matter 
density and therefore they contribute in the same way to the rotational velocities. We also 
define the function 

X(x,v) = -^-[er%) -erf (a;)] (Al) 

where 

erf (x) = — T exp(-t 2 ) dt . (A2) 

y IX JO 

The normalization constant iV of the Maxwellian DF is given in Eq.([5]). 



1. Non rotating model 



N( 3 \V2 

Z{V min, Vq, fesc) - — (t^^ ) X 

X(aw - V, x mm + r]) -2r] exp(-2; 2 ) x min < z - rj 

X(x min -i],z) -exp(-z 2 )(z + r]-x min ) z - rj < x min < z + r] (A3) 
x mill >Z + 7] 
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2. Maximally co rotating Lynden Bell model 

In this case there are two possible situations, depending on the relative magnitude of v Q 

and v esc : 

i) 2 7] < z 



-^(^min; ^0; ^esc 



N 



3 \ 



1/2 



X 



// \'2~ir- 1 



— exp(—z 2 )(r] — z + \?z 2 + 

X min + ??) + exp(r7 2 )x( ^min, V^ 2 + ?? 2 ) 

-exp(-2 2 )(r] - z + V-2 2 + V 2 ) 
exp(r] 2 )x(x min , y/z 2 + r] 2 ) - exp(-z 2 )(^z 2 + rj 2 - x min ) 




< r\ 



V < Xmm < Z-T) 

Z-7] < X min < y /z 2 + rj' 

Zmin > \JZ 2 + T] 2 



(A4) 



ii) 2 > z 



^(^min, ^esc) — 2 
Xv^min Tj^ X 



N 



Gttw 2 ) 



1/2 



7] \2tTW 2 , 

' V) + X(z, 0) + exp(r] 2 )x(v, V z ' 2 + V 2 ) 



min 

— exp(— z 2 )(i] — z + \/ z 2 + rj 2 ) 

min 

— 77, 0) H- exp(?7 2 )x(^, V^ 2 + ^ 2 ) 

- exp(-^ 2 )(v^+V - x min ) 
exp(r] 2 )x(x min , y/z 2 + r] 2 ) - exp(-z 2 )(y/z 2 + rj 2 - x min ) 




•^min ^ Z Tj 
Z-7] < X min < T] 

< y/z 2 + rj 2 

^min > V-2 2 + V 2 



(AS) 



3. Maximally counter rotating Lynden Bell model 



N 



3 \ 



1/2 



// \2nw 2 J 
' x(0, z) + exp(r] 2 )x(Vz 2 + V 2 , v) 



x 



— exp(— z 2 )(r] + z — y/z 2 + ry 2 ) 



X(^min — 77, ^) + exp(?7 2 )x( v / ^ 2 + ?? 2 , x min ) 

-exp(-z 2 )(r] + 2; - y/z^Trj 2 ) 
X(^min -r),z)- exp(-z 2 ) (z + rj- x min ) 




•^min ^ ?! 



< y/z 2 + r] 2 

Vz 2 + V 2 < X mill <Z + 7] 
•^min ^ Z -\~ Tj 



(A6) 
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TABLES 



TABLE I. Threshold energy E^t, quenching factor Q and resolution parameters a and b for 
the detectors considered in the text. 



Detector 


E* (keV) 


Q 


Resolution parameters 


DAMA/Nal 


2 


Na : 0.31 
I : 0.09 


a = 
b = 0.579 


Ge/Neuchatel 


1.5 


0.25 


a = 
6 = 0.17 


Ge/Twin 


4 


0.25 


a = 
& = 0.05 


DAMA/Xe 


13 


0.44 


a = 0.056 
b = 1.191 
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FIGURE CAPTIONS 



FIG. 1. Maxwell-Boltzmann velocity DF in the local frame (arbitrary units). Solid, dashed and 
dotted lines refer respectively to the non-rotating, maximal co-rotating and maximal 
counter-rotating case. 

FIG. 2. Differential rate on a Ge detector as a function of the electron-equivalent energy E ee . 
The solid, dashed and dotted lines refer respectively to the non-rotating, co-rotating 
and counter-rotating case. The upper lines refer to m x = 60 GeV and a = 10~ 9 nbarn 
The lower curves are for m x = 20 GeV and do = 10~ 10 nbarn. The solid (dot-dashed) 
histogram represents the 90% CL. upper limit on the counting rate of the Neuchatel 



(Twin [I4|) experiment. 



FIG. 3. The 90% C.L. upper limit on the quantity £&s™\lr° n ^ as a function of the WIMP mass 
m x for the Ge detectors |TZ|-|r5|]. The solid, dashed and dotted lines refer to a non- 
rotating, co-rotating and counter-rotating halo, respectively. 

FIG. 4. Differential rate on a Nal detector as a function of the electron equivalent energy E, 



CC j 



for a WIMP of m x = 60 GeV and o"o = 10~ 9 nbarn. The solid, dashed and dotted lines 
refer to the non-rotating, co-rotating and counter-rotating case, respectively. The 
histogram represents the 90% C.L. upper limit on the counting rate obtained from the 



DAMA/Nal Collaboration p 



FIG. 5. The 90% C.L. upper limit on the quantity C,c^.^ on ^ as a function of the WIMP mass 
m x for the DAMA/Nal detector [jnj. The solid, dashed and dotted lines refer to a 



non-rotating, co-rotating and counter-rotating halo, respectively. 

FIG. 6. Relative deviation R = [^o'^^ on ^} ro t / [^a^^ onS) ] non ro t between the 90% C.L. upper 
limits in the case of a rotating and a non-rotating halo, for the DAMA/Nal detector. 
The dashed (dotted) line refers to the co-rotating (counter-rotating) halo model. 

FIG. 7. The 90% C.L. upper limit on the quantity £cTg™ 1 f r con ' ) as a function of the WIMP mass 



m x for the DAMA/Xe detector JT8[. The solid, dashed and dotted lines refer to a 



non-rotating, co-rotating and counter-rotating halo, respectively. 
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